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ABSTRACT 

The inverse Compton catastrophe is defined as a dramatic rise in the luminosity of inverse 
Compton scattered photons. It is described by a non-linear loop of radiative processes that 
sets in for high values of the electron compactness and is responsible for the efficient transfer 
of energy from electrons to photons, predominantly through inverse Compton scatterings. We 
search for the conditions that drive a magnetized non-thermal source to the inverse Compton 
catastrophe regime and study its multi-wavelength (MW) photon spectrum. We develop a 
generic analytical framework and use numerical calculations as a backup to the analytical 
predictions. We find that the escaping radiation from a source in the Compton catastrophe 
regime bears some unique features. The MW photon spectrum is a broken power law with a 
break at ~ m e c 2 due to the onset of the Klein-Nishina suppression. The spectral index below 
the break energy depends on the electron and magnetic compactnesses logarithmically, while 
it is independent of the electron power-law index (,v). The maximum radiating power emerges 
typically in the y-ray regime, at energies ~ m e c 2 (~ y m axf«eC 2 ) for s > 2 (s < 2),where y max is 
the maximum Lorentz factor of the injected electron distribution. We apply the principles of 
the inverse Compton catastrophe to blazars and y-ray bursts using the analytical framework 
we developed, and show how these can be used to impose robust constraints on the source 
parameters. 
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1 INTRODUCTION 

The inverse Compton catastrophe is a non-linear loop of processes 
that is closely related to the synchrotron self-Compton (SSC) e mis¬ 
sion when electron cooling is taken into account (e.g. iLongairl 
1201 lh . If a magnetized, non-thermal emitting source, is sufficiently 
compact, then it is possible that the energy density of the radiated 
photons dominates over the magnetic one, which causes the elec¬ 
trons to lose energy mainly by inverse Compton scattering rather 
than synchrotron radiation. In such a case, a runaway process is 
possible: low energy (e.g. radio) photons produced by synchrotron 
radiation are scattered to higher energies (e.g. in X-rays) by the 
same relativistic electron population. As the energy density of these 
high energy photons is larger than that of synchrotron photons, the 
electrons suffer even greater energy losses by up-scattering them to 
even higher energies, e.g. in y-rays. In turn, these have greater en¬ 
ergy density than the X-ray photons, and so on. This runaway pro¬ 
cess leads, on the one hand, to an increase of the electron cooling 
rate and, on the other hand, shifts the radiated photon power from 
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low to high frequencies. In particular, the radiated power is shifted 
to the y-ray regime, up to energies of the order of the maximum 
electron energy. Pushing this idea to its limits, the non-linear loop 
of processes can lead to complete and very fast electron cooling 
and from this, the term inverse Compton catastrophe is coined. 

An important concept in the study of the inverse Compton 
catastrophe are the multiple up-scattered photon generations. One 
can label the photon generations produced in the source according 
to how many scatterings they have undergone after production by 
the synchrotron process. The synchrotron photons are defined as 
the zeroth order photon generation to be produced in the source, 
which, in turn, serve as the seed photons for the first inverse Comp¬ 
ton scattered photon generation; similarly, photons belonging to the 
i-th generation are the seed photons for the i + 1 generation. The 
production of such higher order photon generations can take place 
even if the source is optically thin to Thomson scattering, and in 
this case the dominant loss of photons belonging to a certain gen¬ 
eration (with i > 0) is the escape from the source within a cross ing 
time (see also lBiornsson & Aslaksenll200d : lTsang & Kirkll2007l) . 

The inverse Compton catastrophe has been mainly studied in 
the context of compact radio sources (e.g. lKellermann et al.ll 19691 : 
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lReadheadll994l : lBiornsson & Aslaksenl2000l : lTsang & Kirkl2007t) . 

It has been initially invoked to explain the observed upper limit 
in the distribution of brightn ess t emperatures (7b r ) from powerful 
extragalactic radio sources iKellermann et alJI 1 9691) . Later it was 
suggested that the observed cutoff in the T b r distribution was not 
related to the inverse Compton catastrophe, since the latter sets an 
upper limit of T/ f ~ 10 12 K, which lies above the observed cutoff 
iReadheacjl994f) . Yet, T b c r c remains a theoretical upper limit, against 
which the observed brightness temperatures, e.g. derived from the 
intra-day varia b ility i n radio cores of blazars, are still tested (e.g. 
Fuhrmann et al120081 and references therein). 

The efficient cooling of non-thermal electrons due to inverse 
Compton scatterings has been also studied in a different context, 
i.e. that of unmagnetizecf] sources with high optical depths for 
photon-photon ab sorption and Comp t on scattering. A large volume 
of literature (e.g . iFabian et alj 19861: Lightman & Zdziarski 1987t 


Svenssonl Il987l : lzdziarski 19881 : Zdziarski et all 1 990l : Svenssonl 

1994) deals with the problem of pair (saturated or not) cascades 


that are (i) initiated by absorption of high-energy photons on a soft 
(fixed) photon field and (ii) mediated by efficient Compton cool¬ 
ing of the pairs. The main goal of these studies was to determine 
the properties of the escaping radiation under the modification of 
non-thermal pair cascades and explain, among others, the univer- 


tive galactic nuclei (e.g. Bonometto & Rees 

1971; Rothschild et al. 

198J 

Kazanas 19841; Zdziarski et al. 1990; 

Svensson 1994}). 


Here we consider the production of multiple photon genera¬ 
tions due to inverse Compton scattering in magnetized, optically 
thin to Thomson scattering, sources, and in the regime where pair 
cascades are negligible. The physical setup is si milar to previ¬ 
ous s tudies of the inverse Compton catastrophe (e.g. lTsang & Khj 
I2007I) . yet our main goals differ. These are summarized in the fol¬ 
lowing: (i) study in a systematic way the transition from linear syn¬ 
chrotron to non-linear SSC cooling under the presence of higher 
order up-scattered photon generations; (ii) develop a general an¬ 
alytical framework that is easy to be used and can be applied to 
different types of sources; (iii) determine analytically the condi¬ 
tions, such as the electron compactness, above which the higher 
order SSC photon generations dominate the energetic output of the 
source; (iv) derive an analytical expression for the spectral index 
in the limit where the multi-wavelength photon spectrum is deter¬ 
mined by higher order photon generations; and (v) use numerical 
calculations to back up our analytical predictions. 

This paper is structured as follows. We present our assump¬ 
tions and our analytical treatment in Sect.[2]and continue in Sect. [3] 
with a presentation of our numerical results and a comparison 
against the analytical predictions. In Sect. [4] we present indicative 
applications of the inverse Compton catastrophe in the context of 
high-energy emitting astrophysical sources. We conclude in Sect.[5] 
with a discussion of our results. 


be homogeneous and magnetized, and to contain a population of 
relativistic electrons, which is constantly being replenished. In this 
framework, we will derive (under certain simplifying assumptions) 
analytical expressions for the conditions leading the source in the 
inverse Compton catastrophe limit. For this purpose, we will take 
into account not only synchrotron cooling of electrons, but also 
cooling on synchrotron and higher order inverse Compton-scattered 
photons. In order to keep the analysis as general as possible, all 
of our results will be expressed in terms of three basic quantities: 
the electron, the magnetic and the external photon compactnesses, 
where the compactness is a dimensionless measure of the respec¬ 
tive energy densit)@ and ext ensively used in studies of the emission 
from compact sources (e.g. ZdzjarekL&Yjghtman 1985; Sv enssonl 
ll987l : lLrghtman & Zdziarskill 19871 : Mastichiadis & Kirkll 1995h . 


2.1 Assumptions 

We consider a homogeneous spherical region of radius r b and vol¬ 
ume 14 = 4nd /3 that immersed in a magnetic field of strength B, 
and moves relativistically with Lorentz factor T. As the relativistic 
motion of the emitting region will not become central to our anal¬ 
ysis until much later, the only reference frame of interest, at this 
point, is the rest frame of the emitting region, where all the follow¬ 
ing calculations will be performed. 

For the purposes of the analytical treatment, we assume that 
mono-energetic relativistic electrons with Lorentz factor yo> 3 are 
being injected in the emission region with a constant rate Q 0 , while 
injection of a more realistic power-law distribution will be consid¬ 
ered in the numerical calculations. We further assume that electrons 
may physically escape from the emission region at the (energy in¬ 
dependent) characteristic timescale f e esc = r b /c. This is also equal 
to the typical timescale for adiabatic losses, at least for a spherical 
source expanding close to the speed of light. Electrons lose energy 
through synchrotron radiation and inverse Compton scattering on 
the internally produced synchrotron photons (SSC) and possibly, 
on external photons (external Compton or EC). In order to keep the 
analysis as general as possible, we will not specify the origin of 
the external radiation until much later (Sects. [23l andl4li. The only 
property of the external photon field that is currently important for 
our analysis is its energy density u ex as measured in the rest frame 
of the emitting region. 

We also approximate the emissivity of synchrotron and inverse 
Compton scattering by a d-function centered at (dimensionless) 
photon energies x syn = by g and jc ssc ,i = (4/3)'£y 2(,+1) , respectively. 
Here, x SSCii is the energy of a photon that has been Compton up- 
scattered i-times, b = B/B a with B being the comoving magnetic 
field strength and B cr = 4.4 x 10 13 G. Finally, in the analytical treat¬ 
ment we consid er only optically thin to synchrotron self-absorption 
(ssa) cases (see lTsang & Kirkll2007l . for the role of ssa in Compton 
catastrophe cases). 


2 ANALYTICAL RESULTS 

We adopt a rather abstract framework for the description of the 
source, according to which, the emission region of the source is 
assumed to be embedded in an external (fixed) radiation field, to 


2.2 Steady-state solutions of the electron kinetic equation 


The evolution of the electron distribution n c is described by the 
kinetic equation 


dn e (y, t) n e (y, t ) 

dt fe,esc 


Qeiy , 0 + £*(y. 0, 


(l) 


1 A first analytical treatment of pair cascade in the sync hrotron self- 
Compton model can be found in i Zdziarski & Lightmanl 119851) ; the problem 
of pair cas cades in magn etized sources has been studied also numerically, 
later on, by lCoppil 1 1992b . 


2 The (dimensionless) compactness is traditionally defined as ~ 
(TjL/4nm e c^R , where L and R are the luminosity and size of the system 
iHertericHl 19741) . 
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where Q e and £, e are the injection and energy loss operators for 
electrons, respectively. These are defined as 


Qe 

= QoS(y - yo)H(t), 

(2) 

L e 

4 °tc d r 2 ! 

= a 2 a r*e(y.0wtot , 

3 m e c z oy L J 

(3) 


where H(x) is the Heaviside function and u tot is the total energy 
density. This is given by 


JV T -1 

“tot = ng + U cx + U s yn + ^ I itssc.i ■ (4) 

i=l 

The last term takes into account the energy density of successive 
SSC photon generations, starting from the first one (i = 1) that is 
the result of the synchrotron photon up-scattering (i = 0) and ter¬ 
minating at the (IVt - 1) photon generation. These photons are, in 
turn, the seeds for the last SSC generation to be produced in the 
Thomson regime. SSC photons belonging to the (Vx-th generation 
appear in the electron’s rest frame with energy » m e c 2 , and thus, 
their up-scatter will take place in the Klein-Nishina regime, where 
the cross section for scattering is greatly suppressed. This also ex¬ 
plains the absence of n SSCiNl . from the electron cooling term. The 
highest order of up-scattered photons in the Thomson regime, (V T , 
is given by 


log (1/fryo) 
log(4y2/3) 


(5) 


where the square brackets denote the closest (from below) integer 
value of the enclosed expression (see Appendix [A] for the deriva¬ 
tion). 

In principle, one should also include in the equation of u, ot a 
term related to the repeated inverse Compton scatterings of exter¬ 
nal photons, similar to the last term of the r.h.s. in eq. 0. It can be 
shown, however, that for yo ~ 100 and typical values for the en¬ 
ergy of external photons (e.g. 1 eV), even the first inverse Compton 
scattering of the EC photons takes place deep in the Klein-Nishina 
regime (see Appendix [B}. In this case, one can safely ignore the 
term related to the energy density of EC photons. However, this 
does not apply if the injected electrons are less energetic, e.g. their 
Lorentz factor is ~ 10. Thus, for completeness reasons, we derive 
in Appendix [B] similar relations and constraints to these presented 
in the following paragraphs (Sects. 2.2.1-2.2.2) after including in 
eq. 0 the energy density of the once EC scattered photons (« EC ). 
The interested readers can compare how the inclusion of u EC alters 
the results presented in Sect. [2] 

At this point, it is useful to introduce the compactness of the 
magnetic and photon fields that appear in eq. 0. We define the 
compactness of a field, which is a dimensionless measure of its 
energy density, as 


o- T r h u j 
m e c 2 


( 6 ) 


where the subscript j takes the values 'B’ (magnetic field), ‘ex’ 
(external photon field), ‘syn’ (synchrotron photons) and ‘ssc,i’ (i-th 
generation of SSC photons). Similarly, one can define the injection 
compactness of electrons, which is given by 


^ nj = QtC = rgcrxgoyo 
e Anr^m^c 2 3c 


(V) 


In most astrophysical applications, the total energy density 
n tot is dominated by the first two or three terms, namely by the 
energy density of synchrotron photons and that of the magnetic 


field (or of ambient photon fields). In such cases, a nalytical time - 
depend ent soluti ons of eq. 0 can be found (see lKardashevlll9621 : 
IZacharias & Schlickeisejf20ld , l2012l ). However, if 1 n ssci be¬ 
comes comparabl^fto the sum of the first three terms, then even the 
derivation of steady-state solutions of eq. 0 becomes challenging, 
as the energy density of each SSC generation depends on n e itself. 
In such a case, one can still derive analytical steady-state solutions 
of eq. 0, but only in two limiting regimes, depending on which is 
the dominant term in the kinetic equations. 

In order to define the aforementioned regimes, we introduce 
the ratio £ of the cooling and escape timescales of electrons with 
Lorentz factor yo: 


f = 


4cr T /- b n tot yo 44 t yo ’ 
where 4* is defined accordingly to eq. 0 


We 


distinguish between two regime^ 


( 8 ) 


• If £ > 1, then electron cooling is negligible and the steady- 
state electron distribution is determined by the balance between the 
escape and injection processes. The steady-state electron distribu¬ 
tion is then simply n c (y)/fe,e sc = Qe (y) or equivalently n e (y) = 
k s S(y - yo). Borrowing the terminology from GRBs, we will refer 
to this regime as slow cooling regime (see e.g. lSari et al.ll 1998h . 

• If it, < 1, then electron cooling becomes important. In this case, 
the cooling Lorentz factor of electrons, which is defined by £(y c ) = 
1, is smaller than the injection one (y c < y 0 ) and the steady-state 
electron distribution has the well-known y~ 2 power-law distribution 
(for mono-energetic injection), i.e. n c (y) = k { y~ 2 , for y c < y < 
yo; this corresponds to the so-called fast cooling regime (e.g. I Pi rail] 
120041) . 


The normalization factors k s and kf, whose explicit form is 
presented in the following paragraphs, depend on the injection rate 
Qo and on the energy densities of the various fields. In the para¬ 
graphs that follow, we will determine the slow and fast cooling 
regimes, and present the respective analytical expressions for the 
compactnesses f y 


2.2.1 Solutions in the slow cooling regime 


In the slow cooling regime and for an energy independent escape 
timescale, the steady-state electron distribution has the same en¬ 
ergy dependence as the source function and the normalization k s 
depends only on the injection rate as k s = Qoro/c. Using the defini¬ 
tion of eq. 0 we may write 



OTr b yo 


(9) 


The synchrotron and i-th generation SSC compactnesses are writ¬ 
ten as 


4 y „ = Af B Cyo ( 10 ) 

4c,i = 4(4Cyo)' + \ (11) 

where we used eqs. 0, ( IA7I and jA 101 (see AppendixlAlfor more 
details). The above expressions demonstrate the well-known linear 


3 We will quantify this statement in Sects. 2.2.1-2.2.2 

4 We caution the reader that the slow and fast cooling regimes defined here 

are, strictly speaking, valid only asymptotically, i.e. for £ » 1 and £ 1. 
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and quadratic de pendence of the synchrotr on and first SSC compo¬ 
nents on 4 (e.g. lBloom & MarscheJ 19961) . Since 4c,i oc a cu¬ 
bic relation between the synchrotron and second SSC photon com¬ 
ponents is expected, which might be of interest for flaring events in 
high-energy emitting blazars lAharonian et al.ll2007l) . 

For 4 < l/4yo, expressions eqs. ( Hot and ( fTTl i suggest that 
4c, 4yn < £b , and the total compactness of the source, in this 
case, is Aot ~ + 4- The spectral energy distribution (SED) of 

the source will be synchrotron dominated if 4 > 4 and EC domi¬ 
nated, otherwise. 

The compactness of successive SSC photon generations will 
become progressively higher, namely the source may enter the in¬ 
verse Compton catastrophe regime, only if 4 > l/4y 0 - In this 
case, the following ordering among the compactnesses holds: 


2.2.2 Solutions in the fast cooling regime 


In the fast cooling regime the steady-state equation of electrons is 
written as 


d_ 

dy 


y 2 «e(y) 


mb + w ex + u S s 


Nj —1 

i + ^ ^ Mssc,i 
/=1 


3 ge(y)ro e c 2 

4 CTjC 


( 20 ) 


An Ansatz for the solution n e of the above equation is n c = fc f 7 p - 
By substituting it in the above equation we find that p = 2 and that 
kf satisfies the following algebraic equation 


Nj —1 


/ex + ak { + ^ (ak f ) ,+l 


9C/4 

4cr T r b y 0 ’ 


( 21 ) 


where 


e 


ssc,Nx~l 


» 4„ > 4- 


(12) a = 4cr x r b y 0 /3. 


( 22 ) 


The total compactness is, in a good approximation, given by 

4t ~ 4 + 4 (4C j y 0 ) iVT • (13) 

Depending on the relative ratio of the two terms appearing above, 
the SED is expected to be dominated either by the EC component 
or by the highest order SSC generation. Thus, even for an external 
photon field with high compactness, the production of higher order, 
and progressively more luminous, SSC photon generations cannot 
be suppressed for a sufficiently high electron compactness, namely 
4 > (l/4y 0 )(4/4) 1/A,T ' This will be demonstrated in more detail 
with numerical results in Sect. [3] 

The results we derived so far are valid as long as the assump¬ 
tion of slow cooling is valid, and as such an a posteriori check is 
necessary. The expressions <H-Cl 3 are valid as long as f > 1 or 
equivalently 4t < 3/4yo, where 4t is given by 


4 + 4, 


4 < r- 

e 4yo 


(14) 


4 + 4(44 nj rof T , C>^- 

By combining the above, we find that the slow cooling solutions 
are valid, if the parameters satisfy the following conditions: 


Given that there is no general formula for the roots of a polynomial 
with degree higher than four, which corresponds to Nf = 3, explicit 
solutions of eq. HD can be found in limiting cases only, which we 
list below. 


• Synchrotron and inverse Compton cooling (on synchrotron 
and external photons) 

The main energy loss channels for electrons are synchrotron cool¬ 
ing and inverse Compton scattering on external and synchrotron 
photons, i.e. 4 + 4 + 4yn > 2/I T i * 4c,i- In this case, k { satisfies a 
second degree polynomial equation, which has the solution 


k f 


/ex 

2 a 


-1 + 


[. 12C/4' 

w~ 


(23) 


The above expression for 4 < 4 simplifies into the one obtained 
for SSC electron cooling (see also lPetropoulou et al.l 02013B ). Using 
eqs. © and <A7t . we find the synchrotron photon compactness to 
be 


Cyn — 4«4 — 4 


-1 + 


f, -C.A. 

V “srj 


while the compactness of z-th SSC generation is simply 


(24) 


4 < — and 

4yo 

, _3_ 

/ex < 44yo ’ 

or 

C > -and 

4yo 

where 


/ex = l 



(15) 

(16) 

(17) 

(18) 

(19) 


The quantity / ex characterizes the relative strength of the external 
photon energy density relative to the magnetic field energy density, 
and as such the importance of EC relative to synchrotron and SSC 
emission. It is one of the two parameters that we use to characterize 
the relevant parameter phase space. Relations ([J5} (dll and CD- 
G! can be used to define two distinct regions in the slow cooling 
parameter space. We will refer to them as S 1 and S 2 , respectively. 


4c,i = 41^-j • (25) 

Expressions d and J25b demonstrate how the inclusion of elec¬ 
tron cooling affects the dependence of the photon compactness of 
various components on the injected electron luminosity. For exam¬ 
ple, the synchrotron luminosity does not depend linearly o n 4 J a s 
in the slow cooling regime (e.g. lBloom & Marschenll996h , unless 

124 / 4/4 « 1 . 

We remind the reader that the expression for kf was derived under 
the assumption of akf < 1 (see also eq. Pit ) that allowed us to drop 
the sum over the energy densities of higher-order SSC photons. The 
condition ak { < 1 translates into the following constraint 

4 

/ex >3-^-1. (26) 

t-B 

• Cooling on higher order SSC photon generations 
Here we consider the case where the system is driven in the 
Compton catastrophe regime and the dominant energy density is 
Z4 T i * ^ssc.i ~ 4 c,n t - 1 ■ The energy density of external photons may 
dominate over that of SSC photons, unless 

( 1 ar*f 4 > /ex - 1. 


(27) 
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In this limit, the solution of eq. HD is given by 



a 



(28) 


The above expression is valid as long as the condition ( I27l > holds, 
or equivalently if 


( , »inj x'Vt/A't + D 

/e x<1 + (^j , (29) 

where we substituted expression d28l) into eq. HU. This is the sec¬ 
ond constraint that we derive in the fast cooling regime. 

Using the same equations from the Appendix [A] as before, we 
find that 

( 3 cV /(W t+1) 

4y„ = 4 (3°) 

/ , „inj n(i+1)/(JV t +1) 

4c, i = 41-^-J ■ (31) 


Before closing this paragraph we complete the set of constraints for 
the fast cooling regime. For this, we require the condition £ < 1 or 
4 > 3/4yo to hold. Using the adequate expression for the total 
source compactness, i.e. 


4 + 4 


jnj 

/ex >3^-1 


/ o/>' n j \ At+I / /i* 1 j V 

[ , /ex< 1 + (3%-) 

we find the following constraints: 


, ^ ,3 C , 

/ex > max [ —-,3--1 


^,inj \N T /(Nj + 1) 


( C ) 

/cx < l + (3 t ] 


44yo 4 

^inj n^t/Wt+D 


and 


c 


4/ 3 


N p+ 


3 \44y 0 


(32) 

(33) 

(34) 

(35) 


Similar to what we found for the slow cooling regime, the con¬ 
straints ( 133b and 04 M35t divide the fast cooling parameter regime 
into two sub-regions. We will refer to them as regions F t and F 2 , 
respectively. 

We derived the conditions that result in the dominance of 
the inverse Compton emission over the synchrotron one, both in 
the slow and fast cooling regime. All the conditions were ex¬ 
pressed, so far, in terms of the various compactnesses that describe 
the source. The strength of the inverse Compton emission is tra- 
ditionally expressed also in terms of the Compton Y parameter 
IRybicki & Lightman|[l979l ). For completeness, we present the re¬ 
lation between the Compton Y parameter and the various compact¬ 
nesses in Appendix [Cl The interested reader can use this relation to 
express all the analytical expressions derived in Sect.[2]in terms of 
the Y parameter. 


2.2.3 The / ex - C s parameter space 

The regimes defined in the previous paragraphs are illustrated in 
Fig. D] where we plot / ex against 4 for yo = 100 and two values 
of the magnetic compactness, i.e. 4 = 3 X 1CU 5 (left panel) and 
4 = 3 x 10 -4 (right panel). For a given 4. /ex translates into a value 
of 4, and thus, from the / ex - 4 plot one can read the values of all 


the basic compactnesses that describe the source. For 4 < 1 the 
source is found to be optically thin to Thomson scatterings, where 
our analysis is valid (see also Sect. l3.3l i. 

On both plots, there are some regions marked as S 1,2 and F 1 , 2 , 
which were respectively defined in Sections 2.2.1 and 2.2.2. To 
demonstrate the regions of the parameter space where the inverse 
Compton catastrophe occurs, we used the abbreviation C.C.. The 
dashed and solid black lines correspond to the slow cooling regime 
constraints of eqs. G3G! 1 and GD-tm respectively. The con¬ 
straints obtained in the fast cooling regime, which are given by 
eqs. ( 133b and eqs. j34M35t are shown with solid and dashed red 
lines, respectively. Finally, the transition region between the two 
limiting cases in the fast cooling regime (see Sect. 2.2.2), is not 
tractable by analytical means. 

The main conclusions drawn from Fig. Glare summarized be¬ 
low: 

• for a fixed 4. the system may enter the fast cooling regime 
either for high values of 4 (lower right region) or high external 
photon compactnesses (upper left and middle regions). In the re¬ 
gion F 1 , the main channel of electron energy losses is synchrotron 
radiation and EC scattering. In fact, there is a transition from syn¬ 
chrotron to EC scattering cooling, as one moves towards higher 
values of / ex within the F t region. 

• in the F 2 region, electrons lose energy preferentially by the 
SSC process. Moreover, the energy density of higher order SSC 
photon generations becomes larger, and as such, it dominates the 
total energy density which appears in the electron cooling term. 

• the regions S 1 and S 2 of the slow cooling regime occupy the 
lower left and middle part of the parameter space, where both the 
electron and external photon compactnesses are low/moderate. 

• the inverse Compton catastrophe may occur either in the slow 
cooling or in the fast cooling regimes. Higher values of 4 are re¬ 
quired in the latter. 

• even for 4 » 4 (or /ex » 1) the inverse Compton catas¬ 
trophe cannot be avoided, if the compactness of electrons is suffi¬ 
ciently high (see also eqs. O and GD). 

• an increase of 4 results in a more extended fast cooling 
regime, as expected (see e.g. regions F, in the left and right panels 

of Fig.GK 

• the minimum 4, which is required for the inverse Compton 
catastrophe to occur while in the fast cooling regime, decreases as 

oc (T B llN \ 


2.3 Observed luminosities 


Having investigated the different regimes of the available parameter 
space, we proceed with the calculation of observed quantities, and 
in particular, of the bolometric synchrotron and SSC luminosity. 
For a spherical region moving relativistically with a Doppler factor 
O - r 1 (1 -jgcos 0 O bs)~/ where T is the bulk Lorentz factor and 
8 abs is the angle between the beaming direction and the line of sight 
to the observer, the observed synchrotron and SSC luminosities are 
easily derived by the respective compactnesses as: 


L obs 


p Anr^m^c'TT 

(X x 


i = syn, ssc, 


(36) 


where is given by eqs. G 3 -GD for the slow cooling regime and 
by eqs. P4b -(125l> and |30b-<3 1 b for the fast cooling regime. From 
this point on, all quantities noted with the subscript ‘obs’ will refer 
to the observer’s frame, while quantities with an asterisk (★) will 
be measured in the rest frame of the central engine/galaxy. 

The details about the external source of photons, such as its 
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logle* 



Figure 1. Log-log plot of / ex versus tp for yo = 100. Left and right panels are obtained for 4 - 3x 10 -5 and 3 x 10 -4 , respectively. The constraints defined 
by eqs. d 15t - 4 1 6t and i ll7Ml8) are shown respectively with black dashed and solid lines. The constraint imposed by eq. {33) is plotted with a solid red line, 
while the red dashed line corresponds to the constraint given by eqs. {34M35) . The various sub-regions defined in text are marked on the plot as S 1,2 and F 1 , 2 - 
The gray-hatched region denotes the parameter regime in the fast cooling limit, where our analysis is not valid. 


emission region inside 

emission region outside 

emission region far outside 

external photon region 

external photon region 

external photon region 



Figure 2. Illustration of the physical system under consideration. The emission region (yellow blob) is part of a relativistic jet launched by the central engine, 
and lies at a distance r from it. In the rest frame of the central engine, the external photon field is isotropic and spherical with radius r ex (purple region). 
If r < r ex , the energy density of the external photon field is significantly boosted (oc r 2 «ex,*) because of the relativistic motion of the emission region. The 
opposite holds for cases where r»r a , while the effect of Lorentz boosting is intermediate for r > r ex . 


geometry and spectral energy distribution, were not important for 
the derivation of the results presented in Sect. 2.2. Before stepping 
into the calculation of the EC observed luminosity, it is, however, 
useful to specify the so far arbitrary source of external photons. 
We consider the case of an isotropic (in the galaxy rest frame), 
monochromatic photon field with a photon energy and spher¬ 
ical geometry with a radius r ex (see Fig. 0. Its energy density is 
w ex> * = L cx /4nr^ x c, where L ex is its luminosity. If the emission 
region is embedded in the external photon field, the energy den¬ 
sity and the characteristic photon energy of this field in the co¬ 
moving frame of the emiss ion region appear as r< ex =* T 2 u ex i , and 
5 - x ^ ree X ,*, respectively (Dermer & Sch licke isei]fl993l : |d ermeil 
11995h . If, however, the emission region is located at a distance 
r » r ex and is moving away from it with a Lorentz factor T such 
that /-/r ex » T (see e.g. Fig. [2}, the external energy densit y will b e 
then significantly (by a factor of ~ 1/r 2 ) de-boosted (e.g. lPermetl 
Il995h . In this case, the system can be satisfactorily studied by tak¬ 
ing into account only the internally produced photons. Finally, for 
cases where r/r ex > T the boosting effect on the energy density 
and typical photon energies of the external radiation field will be 
intermediate. 


For the obse rved EC luminosity (L™) we follow the approach 
byD erme 3 fl995b - henceforth D95, since the external photon field 
in the comoving frame of the emission region is not isotropic. For 
the slow cooling regime, where the stationary electron distribution 
is mono-energetic, i.e. n c = k s S(y-yo), we start with eq. (5) of D95 
and derive (see Appendix[P]for more details) the observed flux due 
to EC scattering 


F“ s (e obs , n obs ) = F 0 


| £>(l+^ obs ) j 3/2 | 


gpbsfl +Z) 


V /2 

I <5 (gobs _ Sec) . 


( 37 ) 


where f3 = y/l - 1/r 2 , £ec is defined as 

_ -2) £ex,*yo( 1 /^obs) 

£ec = (1 +z)(l +P) ’ 

where ^ obs = cos 6 obs and F 0 is given by 


(38) 


Fq = 


2) 4 V b cr T c 

4 


Uexkyo- 


(39) 


After integrating eq. 07} over e obs and making the approximations 
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#obs — 1/r, yU„bs ~ 1 - l/2r 2 , 1 + /3 - 2 and T) — Y we find 

I°£ * Y^V^cu^yl (40) 

For the fast cooling regime, where the stationary electron distribu¬ 
tion is a power-law, we can use directly eq. (7) in D95. After mak¬ 
ing the substitutions K —» k f and a = 1/2, where a = (p — l)/2, we 
integrate eq. (7) over energies up to e obs ~ S 2 (l + Fobs/To^ex,*/(1 + 
j8)(l + z), and by making the same approximations as before, we 
derive a similar expression to that of eq. {40}: 

L^xY 6 V b o- T cu, x k f y 0 . (41) 

Thus, Lg^ s depends on the cooling regime only through the factors 
k s and k f . 

For the expressions regarding the slow cooling regime we 
made use of eqs. ©-{ED, {36}, and {40}. whereas for the fast cool¬ 
ing regime we used eqs. {23}-{25}. {30}-{3T}, {36}, and {43. Our 
results are summarized below. 

- Slow cooling regime 


j obs 

-^syn “ 

^r 4 i/ b « B (4Cy 0 ) 

(42) 

r obs 

-^sscj - 

^CV b u B (C4y 0 j +1 

(43) 

r obs 

^EC " 

j^r 6 u b M ex (4Cyo) • 

4r b v ' 

(44) 


- Fast cooling regime 

In Section [2.2.21 we derived explicit expressions for the compact¬ 
ness of various components for two limiting cases: 


• Synchrotron and inverse Compton cooling on synchrotron 
and external photons 


j obs 
^syn 


r obs 
^ssc,i 


j obs 
^EC 


^-Y 4 V b u B f e , 
2 r b 


-1+ A 1 + 


12 17 


f 2 

J ex 


(45) 


3c , 

— 

3c , 

—r 6 v b M ex / ex 
8r b 


-1+ a 1 + 


12C/4 


f 2 

J ex 


-1+ A 1 + 


12 C/C 


f 2 

J ex 


(46) 


(47) 


• Cooling on higher order SSC photon generations (CC limit) 


j obs 
^syn 


j obs 
^ssc,i 


j obs 
^EC 


—r 4 v bWl J—'I 


3r b 


l/(iV T +l) 


E B J 


—r 4 v b M B —— 
3r b B { t B 

3c, 


4 A 


r 6 V b Mex 


3^onj sO'+D/^r+l) 

—) 

(£)' 


1/CJVt+I) 


(48) 

(49) 

(50) 


We remind the reader that the above expressions are strictly valid 
in the slow (regions S 1 , 2 ) and fast (regions Fj, 2 ) cooling regimes, 
and one should use them with caution for parameters falling in the 
transition regime (e.g. hatched regions in Fig.0. 

Thanks to the general formalism we adopted, i.e. all the ob¬ 
served luminosities are expressed in terms of the electron, magnetic 
and external radiation compactnesses, eqs. {ED-©)} can easily be 
applied to a variety of astrophysical sources, such as y-ray emitting 
blazars and GRBs. 


2.4 Photon spectrum in the CC limit 


Besides the observed luminosity of the various components that 
build the multi-wavelength spectrum in the Compton catastrophe 
regime, we are also interested in the spectral shape. As already 
described in Sect. [D the inverse Compton catastrophe refers to 
the case where higher order SSC photon generations are progres¬ 
sively more luminous, i.e. L obs . > L obs . , > ... > L° bs . The 
ordering of the emitting components in terms of their luminosi¬ 
ties, with the synchrotron component being the least luminous, im¬ 
plies that the photon spectrum will be a power-law, defined by the 
peaks of successive components, which will extend up to an energy 
~ 5(yo/10) MeV (as measured in the comoving frame). 

Being equipped with all the necessary analytical expressions, 
we can turn now to the spectral index, fl cc , (approximating the pho¬ 
ton spectrum as a power-law, i.e. L° bs oc e~*c) bn the Compton 
catastrophe limit. Assuming that the bolometric luminosity of each 
component is a good proxy for the peak luminosity, we can estimate 
/3 CC through the approximate relation 


j obs / \ 1 -ficc 

ssc,i+l / £ssc,i+l 1 



(51) 


where e ssci ^ (4/3)'fry 2<I+1) /n e c 2 . Using eqs. {43} and {49} for the 
slow and fast cooling regimes, respectively, we find 


lo g( 4 4 ni 7o) 

log(4y2/3) ’ 

log^f/Vs) 

(W T +l)log(4yjj/3)’ 


slow cooling 


fast cooling 


(52) 


Equation {52} reveals the dependence of the spectral slope on var¬ 
ious physical parameters. For example, in the fast cooling regime, 
a decrease of the magnetic field strength and/or of the size of the 
emission region, results in lower t B values, and thus, tends to make 
the photon spectrum harder (smaller ft cc values). Notice that /3 CC 
is independent of the magnetic compactness in the slow cooling 
regime, while higher C s lead to harder photon spectra in both 
regimes. In principle, there is no lower limit on /3 CC . However, very 
hard photon spectra with /? cc < 0.5 would require unphysically 
large l'T due to the logarithmic dependence on the latter. In any 
case, our results cannot be extended to €° s » 1, since in this pa¬ 
rameter regime other processes, such as photon-photon absorption 
and pair cascades, become important and they cannot be neglected 
anymore. 

The spectral index given by eq. ( 152b is obtained for parameters 
drawn from the sub-regions S 2 and F3 in the / ex - parameter 
space. This is exemplified in Fig© where we adopted yo = 100 and 
t B = 3 x 10~ 5 . In the fast cooling regime, the Compton catastrophe 
photon spectra tend to be harder than these of the slow cooling 
regime. Figure[3demonstrates also that the spectra in the Compton 
catastrophe limit become harder as t' e nj increases. However, because 
of the weak logarithmic dependence on E'"* (see color bar), very 
high values of t'T are required for /3 CC <0.5. Finally, we note that 
a choice of a lower y 0 would not have a strong impact on j3 cc (see 
eq. 123), but it might affect the extend of the sub-regions S 2 and 

f 2 . 


3 NUMERICAL RESULTS 

The results presented in the previous section were derived by solv¬ 
ing the simplest form of the electron kinetic equation (see eq. CD) 
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Figure 3. The / ex - parameter space for = 3 x 10 -5 and yo = 100- 
The color coding denotes the values of the spectral index in the Compton 
catastrophe limit (0 CC ). 


and by using a series of approximations. In this section, we present 
the results of a numerical treatment of the problem, which does not 
suffer from the simplifying assumptions used in Sect. [2] albeit its 
own limitations. Our aim is to test whether or not these simplifica¬ 
tions have a radical effect on the results regarding the transition to 
the Compton catastrophe regime, as well as the multi-wavelength 
spectral properties in this limit. 

The numerical treatment of the problem allows us to augment 
the relativistic electron kinetic equation with more processes, such 
as pair injection due to photon-photon absorption, and to make use 
of the full expressions for the synchrotron and Compton emissivi- 
ties in the relativistic limit. In addition, we write an accompanying 
equation for photons, which is coupled to the one for electrons. The 
physical system is then described by 

dn e n s 

tc.csc 

dn y n 7 

ty.csc 

The various terms appearing above describe the following physi¬ 
cal processes: (i) Photon-photon pair production, which acts as a 
source of new particles (<2e 7 ) an d a sink term for photons (_£)7); 
(ii) synchrotron radiation, which acts as an energy loss term for 
electrons (_ C7 ") and a source term for photons (<3) yn ); (iii) syn¬ 
chrotron self-absorption, which acts as a sink term for photons 
(-£y Sa ); and (iv) inverse Compton scattering, which acts as an en¬ 
ergy loss term for electron^ (Xe s ) and an energy injection term for 
photons (<2^' s ) The escape timescale for electrons f e ,esc is set equal 
to the crossing time of the source, i.e. t CT = r\,/c. Note that the 
expressions for the synchrotron and Compton e missivities in the 
relativ istic limit (c.f. eqs. (6.33) and (2.48 ) in lRvbicki & Lightmanl 
( fl979b and iBlumenthal & Gould (l970al) . respectively) have been 
used in eq. ns. For more details about the functional forms of the 
various rates, see Appendix lEl 

Synchrotron self-absorption is taken into account only as a 
sink for photons, while there should be an additional source (heat¬ 
ing) term for the electrons. Although this is not, strictly speaking. 


= LT + XT + or + a 

(53) 


(54) 


5 Both external and internal photons serve as targets for the inverse Comp¬ 
ton scattering process. 


an energy conserving scheme, in all cases we present the absorbed 
energy by synchrotron self-absorption is small, and electron heat¬ 
ing can be neglected. Equations and could not be ap¬ 
plied in the ir pres ent form in the “synchrotron boiler” case, where 
t B » C dGhisellini et alJll988l) . 

As the numerical code makes no distinction between electrons 
and positrons (both determined by eq. ED>. the tre atment of th e 
pair annihilation in the code is approximate (see e.g. lCoppilll992h . 
As long as we consider optically thin cases, i.e. <1-3 (or, 
equivalently r e < 1 where r e is the Thomson optical depth and is 
defined in Appendix C) the pair yield is very low and the annihila¬ 
tion can safely be neglected. The annihilation line starts be coming 
important for 4^ > 10 (e.g. iLightman & Zdziarskill 19871) . where 
our numerical approach becomes invalidated. 

In summary, t he numerical code used in this paper and pre¬ 
sented originally in lMastichiadis & Kirkl d 19951 ) is adequate for the 
description of system containing relativistic electrons (y > 3) with 
compactnesses up to ~ 3. Exploration of a system containing elec¬ 
trons with/3y < 1 and/or with very high c ompactness re quires more 
sophisticated numerical treatment (e.g. Comn|_ l_992l : Istem et aD 
ll995l : lBelmont et alj2008l : IVurm & Poutanenl20091) : this, however, 
lies out the scope of the present study. 


3.1 Comparison with analytical results 

As a first step, we compare the analytically derived values for the 
observed luminosities (see Sect. [23l against the values we obtained 
by numerically solving eqs. (53) and (54). without including the 
terms of synchrotron self-absorption and photon-photon absorp¬ 
tion. In this way, we ensure that any differences found between 
the numerical and analytical results will have to be related with the 
simplifying approximations of the analytical treatment and not with 
the physical processes themselves. 

For the comparison we adopted the fiducial parameter set of 
yo = 100, C = 10~ 2 , B = 10 G, r b = 10 16 cm, T = 10. Substitution 
of the y 0 and B values in eq. (5} results in Nj = 2. The magnetic 
compactness is found to be ( B = 3 x 10 -2 > To test the effect 
of the external photon field on electron cooling and on the multi¬ 
wavelength photon spectra we used three indicative values for its 
compactness, namely t ex = 0,0.1 and 0.4, which correspond to 
equally spaced values of log / ex . In all cases but 4x = 0, the basic 
compactnesses are ordered as £' e ni < C B < 4x- F° r the characteristic 
energy of external photons, we used the fiducial value e ex * = 1 eV. 

Using the analytical expressions of eqs. and 01- 

03 we construct the respective / ex - parameter space, which 
is illustrated in Fig. [4] The parameter values adopted for the com¬ 
parison are shown as red crosses. For all three values of 4x> we 
analytically predict that the system should lie in the fast cooling 
regime (sub-region F,). We note that even for f ex = 0, the magnetic 
compactness {t B x B 2 r b) is high enough as to make the injected 
electrons cool efficiently. This is indicated by the absence of the 
slow cooling regions in Fig. 0 In what follows, we will compare 
the aforementioned predictions with the respective numerical re¬ 
sults. 

First, we verified that the steady-state electron distribution is 
described by n e (y) oc y~ 2 as a result of efficient cooling, for all 
values of C ex - The multi-wavelength photon spectra calculated nu¬ 
merically for the three values of are shown in Fig. [5] Photon 
spectra obtained for 4x = 0,0.1 and 0.4 are plotted with black, red 
and blue lines, respectively. The synchrotron and inverse Comp¬ 
ton components of the total emission (thick line) are plotted with 
dashed and dotted lines, respectively. For illustration reasons, the 
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Figure 4. The / ex - plane for yo = 100, 5 = 10 G and /*b = 10 16 cm 
(£b = 3 x 10 -2 ). The parameter values used in the numerical runs of Sect. [3] 
are shown as symbols: red crosses for the comparison between the analyt¬ 
ical and numerical treatments (for more details, see text) and black circles 
for the transition of the system to the Compton catastrophe regime. The 
red lines and the grey-colored hatched region have the same meaning as 
in Fig. [T] The sub-regions of the fast cooling regime, F\ and Co. are also 
marked on the plot. 



log V (Hz) 

Figure 5. Steady-state multi-wavelength photon spectra calculated numer¬ 
ically for three values of C ex marked on the plot. The total emission is 
shown with thick solid lines, while the synchrotron and inverse Compton 
(SSC+EC) components are shown with thin dashed and dotted lines, re¬ 
spectively. Other parameters used are yo = 100, B = 10 G, ty, = 10 16 cm, 
Cb = 3x HP 2 , = Hr 2 , T = 10 , e ex ,* = 1 eV. The flux values correspond 
to a source with a fiducial redshift z = 0.2. 


external radiation field is not shown. The flux values marked on the 
right vertical axis correspond to a source placed at redshift z = 0.2. 

In agreement with our analysis, we find two SSC photon gen¬ 
erations being produced in the Thomson regime. This is evident in 
the case of pure SSC emission (black lines). The increase of l tx 
has two effects on the shape of the SED. On the one hand, the syn¬ 
chrotron luminosity decreases, while, on the other hand, there is 
a transfer of power from the SSC to the EC component (see e.g. 
red and blue lines in Fig.0. This reflects the transition from syn¬ 
chrotron cooling to cooling due to EC scattering, which is in agree¬ 
ment with the analytical predictions (see Sect. l2.2.3l l. 

Apart from a qualitative comparison, we can compare the 
bolometric synchrotron, SSC and EC luminosities as derived nu¬ 
merically and analytically from the appropriate expressions in 
Sect. [2] namely using eqs. J45M471. The results are listed in Ta- 
bleQ] In all cases, the analytically derived values do not differ more 
than a factor of ~ 2.5 from those derived numerically. Given the 
simplifying assumptions we used to derive the analytical expres¬ 
sions in the first place, we argue that the results are in good agree¬ 
ment to each other. 



Figure 6. Steady-state multi-wavelength spectra obtained for different val¬ 
ues of starting from 10 2 (black lines) and increasing up to III 0 2 (grey 
lines) with logarithmic increments of 0.3. The multi-wavelength spectra ob¬ 
tained for f ex = 0 are over-plotted with dashed lines, for comparison rea¬ 
sons. All other parameters are the same as in Fig. [5] 


3.2 Transition to the inverse Compton catastrophe 

In this paragraph we demonstrate the modification of the multi¬ 
wavelength photon spectra caused by a gradual increase of the elec¬ 
tron injection compactness, which eventually leads the source to 
the Compton catastrophe regime. Same as in previous paragraph, 
we neglect the effects of additional processes, such as synchrotron 
self-absorption, which will be discussed separately in a following 
paragraph. 

We use as a starting point, the fiducial parameter set of the pre¬ 
vious section with £ ex = 0.1 (see Fig. 0 In this case, the external 
photon energy density is the dominant term in the electron cooling. 
This is also reflected to the SED shown in Fig. [5] where most of the 
bolometric photon luminosity is carried by the EC component. We 
will show next that even in this case, higher order SSC photon gen¬ 
erations may have a dominant contribution to the SED and electron 
cooling, for sufficiently high as already pointed out in Sect. [2] 


Since this particular example falls in the fast cooling regime, we 
can use the expressions for T SSCiNt and Tec given by eqs. ( 1491 ) and 
respectively, in order to predict the necessary electron com¬ 
pactness for the highest SSC component to be the most luminous. 
By demanding T ssc ,2 > Tec. we find that t™ 1 > 0.2, yet t'" s < 1. 

We initiate the numerical calculations with = 1CT 2 and 
keep increasing by a factor of two over its previous value until 
we reach = ICC 0 ' 2 . The predicted value for the transition to the 
inverse Compton catastrophe is slightly below the maximum value 
used in the simulation (-0.7 in logarithm). The successive values 
are shown as black circles in the parameter space of Fig. [4] 

The successive photon spectra are shown in Fig. [6] with solid 
lines, while spectra obtained for the pure SSC case, i.e. with 
4 X = 0, are shown with dashed lines for comparison reasons. For 
the particular choice of parameters (see also Fig. 0, we find that 
the EC emission dominates in the y-rays, while the SSC emis¬ 
sion is being suppressed, for low enough values of the electron 
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Table 1 . Observed bolometric luminosities of the synchrotron, SSC and external Compton components of the SED shown in Fig. [5] 



^ex 

= 0 

P 

v e\ 

= 0.1 

P 

v e\ 

= 0.4 


A (a > 

N (b) 

A 

N 

A 

N 

(erg/s) 

8.8 x 10 47 

4.0 x 10 47 

3.2 x 10 47 

2.5 x 10 47 

1.1 x 10 47 

1.2 x 10 47 

C, 1 (erg/s) 

5.1 x 10 47 

2 x 10 47 

7 x 10 46 

7.9 x 10 46 

7.1 x 10 45 

1.5 x 10 46 

L Jc, 2 < er g/ S ) 

3.0 x 10 47 

6.3 x 10 46 

1.5 x 10 46 

hidden (c) 

4.9 X 10 44 

hidden 

L° bs (erg/s) 

- 

- 

2.5 x 10 47 

4.7 x 10 47 

3.2 x 10 47 

7.9 x 10 47 


a Calculated using the analytical expressions given by eqs. E3-E3, which are valid in the fast cooling 
regime (F i). 

b Calculated numerically - see also Fig. [5] 

c The numerical code cannot distinguish between photons with EC and SSC origin, if they are produced 
at the same energy range. For example, if L“ bs 2 < l-'yf ■ the 2nd SSC photon generation will be hidden 
by the EC one. Thus, we cannot numerically calculate its luminosity. 


compactness, i.e. t'P = 0.01 - 0.08. As the electron compact¬ 
ness progressively increases, however, the contribution of the EC 
component to the high-energy part of the spectrum becomes less 
significant. We find, in particular, that it is the second SSC emis¬ 
sion that actually dominates the multi-wavelength emission for 
4T J > 1CT 0 ' 8 (second spectrum from the top). Note that this value is 
in good agreement with the analytical estimate derived above. Fi¬ 
nally, the spectrum obtained for = 10~°' 2 is a typical example 
of emission in the Compton catastrophe regime (see also Fig. 2 in 
iBiornsson & Aslaksenll200ol ) and bears several ‘generic’ features: 

(i) The bulk of the total luminosity is emitted at the high- 
energy part of the spectrum. For the pure SSC case shown in 
Fig. [6] (dashed lines), the SED peaks at approximately e ssc ,N T = 
(r/(l + z))(4/3) iVT &yo A ' T+2 meC 2 - 2 MeV or ~ 5 x 10 2 ° Hz. As al¬ 
ready mentioned in Sect. 2.2, the succession of higher than (V T SSC 
photon generations ceases due to the Klein-Nishina effects. The 
highest energy of the inverse Compton scattered photons should 
be e KN ~ ry 0 m e c 2 For the values numerical example of Fig. [6] we 
find eKN = 0.5 GeV or ~ 10 23 Hz. Notice the abrupt cutoff of the 
photon spectra at this frequency. 

(ii) The SED shows smooth spectral breaks at the transition be¬ 
tween successive emission components. This effect is evident, if 
the electron distribution is mono-energetic or spans over a narrow 
range of energies. However, for a power-law electron distribution 
with index x > 1, the spectral breaks are smoothened out (see next 
section). 

(iii) For even higher compactnesses, namely 0.2 < < 3, we 

verified that the spectrum tends to be a power-law with spectral 
index J3 CC ~ 0.8. This is in agreement with the analytical estimate of 
eq. 113, but it will become more clear with the numerical examples 
that follow (Sect. 3.2.1). 

We have shown that there is a general agreement between 
the analytical and numerical results regarding the electron cooling 
regimes, the luminosities of the various components, the transition 
of the system to the Compton catastrophe regime, and the respec¬ 
tive spectral shape in the Compton catastrophe limit. There are, 
however, several factors that might affect the spectral index / 3 CC , 
such as the distribution of electrons at injection (mono-energetic 
versus power-law) and photon-photon absorption, which redis¬ 
tributes the energy from the high-energy part of the spectrum to 


lower energies through an electromagnetic cascade. As the analyti¬ 
cal treatment has its limitations, in what follows, we investigate the 
role of these factors numerically for a pure SSC case (i ex = 0). 


3.2.1 Effects of the electron power-law index 

So far, we have presented multi-wavelength spectra obtained for 
mono-energetic injection of electrons. However, the spectral index 
/3 CC may depend on the power-law shape of the electron distribution 
at injection. To test this hypothesis, we consider here the case of a 
power-law injection of electrons (with y^n < y < y m ax) i.e. Q t oc 
y- s H(y - y mi „)//(y max - y), with y min = 1, y max = 10 22 . We use 
three indicative values for the power-law index (s = 3,2,1) that 
correspond to electron energy spectra ranging from soft to hard. For 
each case, we calculated the spectral for five values of the injection 
compactness starting with tP = 0.01, and increasing fj." 1 over its 
previous value by +0.5 in logarithm. Other parameters used are: 
B = 10 G, r b = 10 16 cm, leading to 1 B = 3 x ICC 2 . In order to 
isolate the effects that a different power-law index has on the photon 
spectra, we artificially switched off the terms related to photon- 
photon absorption and synchrotron self-absorption from eqs. 153} 
and E3. The results for s = 3,2 and 1 are presented in the top, 
middle and bottom panels of Fig. [7] In all cases, we show the total 
(synchrotron and inverse Compton) emission in terms of photon 
compactness. Different types of lines correspond to different 
marked on the plot. In all cases the steady-state electron distribution 
has cooled efficiently down to y c ~ y min <K y 0 . 

For low enou gh electron compactnesses, w e recover the well 
known result (e.g. iBlumenthal & Gould 1 1970a) that the spectral 
shape of the SED depends on the power-law index of the elec¬ 
tron distribution (see e.g. spectra for = 0.01 in all panels of 
Fig. 0- However, as the system is pushed in the Compton catas¬ 
trophe regime (tf* > 0.3), the SEDs obtain a broken power-law 
shape with a break at ~ m c c 2 as measured in the comoving frame. 
The spectral index below the break energy is universal, in the sense 
that, it is no more related to the power-law index of electrons at 
injection, but depends only on and y 0 (for a discussion on 

the charact eristic Lorentz factor yp in th e case of power-law in¬ 
jection, see IBiornsson & Aslaksenl i2000i) ). However, the photon 
spectra for (comoving energies) > m e c 2 are steep and their spectral 
index is related to the power-law index of electrons at injection. The 
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s=2 



s=l 



Figure 7. Steady-state multi-wavelength photon spectra obtained for a 
power-law injection of electrons with index s = 3 (top), s = 2 (middle) 
and s = 1 (bottom). Different types of lines correspond to the values of 
tP marked on the plot. Also marked is the value of the spectral index J3 CC , 
whenever this can be defined. In all panels, x is the photon energy in m e c 2 
units and £{x) is the differential photon compactness. 

broken power-law shape and the steepening of the spectrum above 
the break energy are both res ults of the Klein-Nishina c utoff effect 
(for a relevant discussion, see lZdziarski & Lambl dl986h ). We note 
that the spectra above the high-energy break should appear even 
steeper, had the photon-photon absorption be taken into account, 
while fi cc would be marginally affected by the development of pair 
cascades (see Sect. [33). 

The numerically derived values of /? cc , whenever they can be 


Figure 8. Same as in bottom panel of Fig.[7]but with the inclusion of syn¬ 
chrotron self-absorption as a sink term for photons (dashed lines). The spec¬ 
tra obtained without synchrotron self-absorption are plotted with solid lines. 

defined, are marked on the plot, and are in relative agreement with 
the value derived by eq. 452}. Moreover, we find that the photon 
spectra in the Compton catastrophe regime become harder as the 
electron compactness increases. In particular, a one order of mag¬ 
nitude increase in leads to a larger /J cc by ~0.12 units. The weak 
dependence of /3 CC on 4” J is in full agreement with the analytical 
prediction (see eq. d52H. On the contrary, a one order of magnitude 
change in fe nJ > before the source enters the Compton catastrophe 
regime, results in significant changes of the photon spectral shape 
(see Fig. [7}. The sensitivity of the multi-wavelength photon spectra 
on before and after the transition to the Compton catastrophe 
regime, is also reflected by the different density of curves in the 
panels of Fig.[7]as f"' 1 increases. 

It is noteworthy that the photon spectra we obtain in the 
Compton catastrophe reg ime are not as hard as those predicted by 
Izdziarski & Lama d 1986l ~). which have a spectral index ~ 0.3 for 
even lower values of the electron compactness. After excluding a 
numerical error as the reason of this discrepancy (see Appendix|F}, 
we argue that the hardness of the photon spectra in the Compton 
catastrophe regime depends also on the nature of the seed photons 
used in the scatterings. Here, we study the effects of repeated in¬ 
verse Compton scatte rings of the synchrotron p hotons by relativis¬ 
tic electrons, whereas IZdziarski & Lambl 119861) considered a fixed 
black-body photon field. 

3.3 Effects of additional processes 

We turn now to discuss additional processes that were not included 
so far and can influence the system. 

(i) Synchrotron self-absorption: for the parameter regime we are 
interested in, this process acts as a loss term for low-energy pho¬ 
tons (see eq. d54b ) and reduces the number density of synchrotron 
photons, which are the targets for the first SSC photon generation. 
It may, therefore, affect the spectral slope in the Compton catas¬ 
trophe regime. In order to test this, we calculated numerically the 
steady-state photon spectra for the same cases shown in Fig. [7} the 
case of s = 1 is exemplified in Fig. [8] In general, inclusion of syn¬ 
chrotron self-absorption makes the multiple SSC components more 
pronounced for low £' e nj , while in the Compton catastrophe regime 
the spectrum is characterized by smaller j8 cc . Qualitatively speak¬ 
ing, this is expected, since synchrotron self-absorption removes 
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Figure 9. Same as in Fig.[7]but with the inclusion of photon-photon absorp¬ 
tion (dashed lines). The spectra obtained without photon-photon absorption 
are plotted with solid lines. For each value of s, the lower and upper curves 
are obtained for tP = 0.3 and 1, respectively. 


low-energy photons that would be, otherwise, multiply up-scattered 
and contribute to the high-energy emission of the source. Quanti¬ 
tatively though, the effect is not significant, since in all cases /? cc 
changed no more by a factor of 0.1-0.2. For example, for s = 1 
and p 1 = 1 (Fig.[^ we obtain y3 cc = 0.8 instead of 0.9 after the 
inclusion of synchrotron self-absorption. 

(ii) Photon-photon absorption: inclusion of this process, makes 
the high energy part of the spectrum, i.e. at x > 1, steeper. This 
is exemplified in Fig. [9] where we show the photon spectra for 
s = 1,2 and 3 calculated with (dashed lines) and without (solid 
lines) photon-photon absorption. For each value of s, we plot the 
spectra for £“ J = 0.3 and 1 are (from bottom to top). Only for $ = 1 
and fe" J = 1 starts photon-photon absorption to have a significant 
impact on the high-energy part of the spectrum. Note, however, that 
the injected secondary pairs do not modify the spectrum below the 
break. Thi s becomes important only for > 10 (see e.g. Fig. 2 in 
ISvenssonl 19941) . In short, for the parameter regime where our anal¬ 
ysis is valid the effects of photon-photon absorption are minimal. 

Finally, other processes that might be important in opti¬ 
cally thick plasmas such as photon down-scattering on elec¬ 
trons, electron-positron annihilation, Coulomb collisions and 
bremsstrahlung turn out to be negligible as the Thomson optical 
depth in all our cases is always much less than unity. 


4 RELEVANCE TO ASTROPHYSICAL SOURCES 

We have shown that the multi-wavelength photon spectra in the 
Compton catastrophe regime obtain a universal power-law shape 
that peaks at observed energies » m e c 2 , where most of the power 
is emitted. In contrast to the photon spectra obtained for low val¬ 
ues of the electron compactness, where one can link the observed 
spectral index with the power-law index of the electron distribu¬ 
tion, in the Compton catastrophe regime this connection can no 
more be established. The same applies for the EC emission, which 
for high enough is hidden below the higher order SSC compo¬ 
nents. These properties make the photon emission in the Compton 
catastrophe unique. 

For the derivation of the analytical relations presented in 
Sect. [2] we adopted the most general framework as possible, i.e. 


we did not attempt to specify the astrophysical source. The results 
regarding the observed bolometric luminosities of the various spec¬ 
tral components are obtained for an emitting region that is relativis- 
tically moving towards the observer and is embedded in a spheri¬ 
cal and isotropic (in the rest frame of the galaxy) external photon 
field. These analytical relations are valid as long as the relations 
w ex — F 2 M ex ,* and «ex — Tfex,* hold. In this section, we apply the 
analytical results of Sect. [2] to the high-energy emitting region of 
either a blazar jet or a GRB jet. 


4.1 Blazar emission 


Assuming that the blazar multi-wavelength emission (from the in¬ 
frared (IR) wavelengths up to y-ray energies) arises from one re¬ 
gion, then the observed spectrum cannot be reconciled with that 
expected in the inverse Compton catastrophe limit. One exam¬ 
ple is the absence of an excessively high X-ray flux from radio- 
loud blazars, whose radio cores exhibit extreme apparent bright¬ 
ness temperatures. This may exclude the possibility of a source be¬ 
ing in the inverse Compt on catastrophe limit (e.g. lOstorero et alj 
120061 : lAgudo et all 1 20061) . If, however, the one-zone assumption 
proves to be wrong, i.e. the IR-optical and y-ray emission are pro¬ 
duced in different locations, then it is still possible to explain the 
high-energy emission in term s of a non-thermal pair cascade (e.g. 
IStern & Poutanenll2006ll2008l) . 

In the following, we adopt the one-zone framework for the 
blazar emission, and, as an indicative example, we focus on blazars 
whose y-ray emission is dominated by the EC component. Thus, 
what follows may be of relevance for flat spectrum radio quasars 
(FSRQs) and low-frequency peaked BL Lacs (e.g. 3C 273 and 
BL Lacertae). In the fiducial example of Fig.[6]we showed that, for 
sufficiently high tP, the emission from the highest order SSC pho¬ 
ton generation may dominate that of EC component. Thus, in the 
case of EC dominated blazar emission, it is interesting to search for 
the conditions that suppress the emission from highest order SSC 
photon generation, and derive information about the properties, e.g. 
size and luminosity, of the external radiation field. 

We present the results for the fast cooling case, since it is the 
regime that ensures high radiative efficiency. The luminosity ratio 
between the high- and low-energy components of the blazar SED 
is usually referred as the Compton dominance q of the source. In 
cases where the EC component dominates the emission in high- 
energies (X -rays u p to y-rays) the ratio is defined as q =s - 

Ucx/ub (e.g. lSikora et al.ll2009l) or q = CtxpB- The condition L^ s 3= 
L° h * results in 

ssc,JNt 


4 q 


> 



Nj/(Nj+l) 


(55) 


where we used eqs. J49t and {30}. Using the constraint of eq. {35} 
we find that 


3C] NtKNt+1) > 3 

£b ) 4f B yo 


(56) 


which when combined with the relation {55} results in r ex < r max 
with r max given by 

Fnax = 3.2 x 10 17 cm (r 2 i 7 _ 1 L dj46 ) 1/2 (r M6 y 0i2 ) 1/2 , (57) 

where T is the Lorentz factor of the emission region. In the above, 
we introduced the notation F x = F/10 X in csg units, unless stated 
otherwise. We also assumed that the source of the external radia¬ 
tion is the reflected blazar’s disk luminosity (L d ) at the broad line 
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region (BLR). The latter, for the purposes of this study, was taken 
to be spherical with radius r ex (see also Fig. □ In the derivation 
of d57b we, thus, assumed L cx = rjL u and used rj^=_ O.lp-i and 
= L a 4 f,/ 10 46 erg/s a s typical values (see e.g. lBottcher et alj2013l : 
INalewaiko et alj|2014l) . 

Relation d57b shows that the inverse Compton catastrophe may 
be avoided, if the external source of radiation, the BLR in our case, 
has a sub-pc scale size (r ex ~ 0.1 pel For a given L cx , a pc-scale 
BLR on the other hand, is dilute enough as to favour the SSC emis¬ 
sion and even lead the source to the Compton catastrophe limit for 
sufficiently high f^ s . In general, the upper limit we derived for r ex 
becomes even more stringent for lower reflection efficiencies 77 and 
less luminous disks. Although a larger emission region and/or the 
injection of electrons with higher Lorentz factors, push the upper 
limit close to the pc-scale, the dependence of r cx on both parame¬ 
ters is relatively weak (r max oc (yor b ) 1/2 ). 

The aforementioned results suggest that the inverse Comp¬ 
ton catastrophe does not typically occur in these sources. Before 
closing this paragraph we note that higher order SSC photon gen¬ 
erations (At > 2 ) may become relevant for y-ray emitting radio 
galaxies, such as Cen A and M 87, even if the sources are not in 
the extreme re gime of the Compton cat astrophe. In the particular 
case of Cen A. IPetropoulou et ail d2014l ~) studied in detail the role 
of the second SSC photon generation on the high-energy emission 
detected from the core of Cen A. 


in the prompt GRB emission. We assume that the GRB emission 
takes place i n the optical l y thin to Thomson scatterings part of th e 
GRB jet (e.g. |Piranll999l : lMeszaros & Reesll2000l : lGianniosll2012h . 
namely 


7"em 


> r T = 10 14 cm 



(58) 


where r em is measured in the galaxy rest frame. The size of the 
emission region, as measured in the respective rest frame, is then 
parameterized as r b - r em /r > lO^Li^T ) 4 cm. An estimate of the 
magnetic field in the emission region is given by 



(59) 


where eg denotes the ratio of the Poynting luminosity to the jet ki¬ 
netic luminosity. The magnetic compactness of the emission region 
is written as 


t^B ~ 


O' TfBLk 

87™ c c 3 r em r 3 


(60) 


The external radiation field, which will be up-scattered by the elec¬ 
trons in the emission region of the GRB jet, is left undefined, 
apart from the fact that is assumed to be spherically symmetric and 
isotropic in the rest frame of the explosion. We also assume that the 
emission region is embedded in it: 


4.2 GRB emission 

Amongst the various proposed mechanisms for the GeV emis¬ 
sion of GRB si , the inverse Compton scattering has been dis- 


'em ” 10 Cm 5; A, 


Lk,53 

r 3 

1 i 


(61) 


This can be considered as the most efficient case of a more gen¬ 
eral scenario where r cm > r ex , since the energy density of exter- 


cussed in several studies (e.g. 

Dermer et al. 20001; Guetta & Granot 

In the limit where r em » r ex , the energy density of external pho- 
tons will appear in the comoving frame of the ejecta de-boosted 

200l;|Granot & Guetta 20o| 

Pe'er & Waxman 20oJ Beloborodov 

2004 Fan & Piran! 200"d Gupta & Zhang) 20071; Ando et al. 20081; (Dermer & Schlickeiser 1993J). 

Fan et al. 2001$). External inverse Compton from various photon 

Two are the observational constraints that we will use in our 


sources has been also considered by man y authors (e.g. [Shavivj analysis: 
Il995l: IShaviv & Darll 1 995l: iGiannio i l2008 ) beginning with lShemil 
( 19941) , while some of its implications when applied as an up- 


sc attering process of the hot cocoon radiation hav e been addressed 
bv lToma et alj d2009h and lKumar & Smootl d20 1 4l) . 

Regarding t he (sub)MeV GRB emission, the volume of litera¬ 
ture is huge ( e.g. |Piranll2004| Pe'er et al fcooel t lBeToborodovlElIol : 


IVurm et al. boi lH G ianniosll 2012 l) . with each study reproducing at 

a certain degree the basic observed properties. Here, we test the al¬ 
ternative hypothesis of an external Compton origin, since it offers a 
physical setup for illustrating the principles of the inverse Compton 
catastrophe and the generality of the analytical framework we de¬ 
veloped in Sect. [2] In what follows, we will show how the inverse 
Compton catastrophe alone, can be used to reject this hypothesis. 

We begin with a short description of the GRB jet model that 
we will use to accordingly define quantities, such C B and r b , that 
appear in the analytical expressions of Sect. [2] Let us consider a 
GRB flow of kinetic (isotropic equivalent) luminosity Z* and bulk 
Lorentz factor T. When the jet reaches a distance r em a substan¬ 
tial fraction of its luminosity is dissipated internally leading to ac¬ 
celeration of electrons. A fraction of the dissipated energy results 


6 Interestingly, the B LR radius is typically found to be 
INalewaiko et al J 2 OI 4 L and references therein). 


0.1 pc (see 


(i) The pea k ene rgy of the GRB spectrum £ P)0 bs ~ 0-2 MeV 
(Gru ber et alj 12014) .which can be used to express the energy of 
the external photons as 


10 eV 


1 e p.°t® ) | 


| 1 

10.2 MeV T 

l 2 J 

'r 2 rg 


(62) 


Because the energy of external photons appears boosted in the co¬ 
moving frame by a factor of T, the above relation sets a strong con¬ 
straint on the product y^e eXti ,, at least in the regime of yo > 3 where 
our analysis is applicable. 

(ii) The GRB energy flux in the energy range 10 keV-1 MeV is 
F° bs ~ 3 x 10 ~ 7 erg cnf 2 s~* liOruber et al.ll2014t) . 

In the scenario of an EC origin of the prompt GRB emission, we 
thus require L = 47 rD[F° bs = 1.5 x 10 51 erg/s, for D L = 6.7 Gpc 
(z = 1). In what follows, we normalize Lg bs with respect to 
10 51 erg/s. 

These observational constraints in addition to the principles of the 
inverse Compton catastrophe can be used to exclude an EC origin 
of the (sub)MeV GRB emission. We demonstrate this in detail for 
the slow cooling regime even though one expects the emitting re¬ 
gion to be fast cooling. The expressions are less complicated in this 
regime and similar conclusions can be reached for the fast cool¬ 
ing regime, although a lengthier derivation is requirecfl Solving 


Fermi-LAT detections dAckermann et al. 

2013b. it has been already ob- 

served with CGRO/EGRET (Hurley et al. 

1994S; 1 Gonzalez et al.l2003h. 8 In the fast cooling regime, the inverse Compton origin of the ~MeV GRB 
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eq. ( 162b with respect to y 0 and substitution into eq. ( 1441) results in 


10 43 erg/s 



where we made also use of eq. ( 161b . 

Relations ED, {62} and ( 163b describe an external source of 
photons with the following properties: r ex > 1500R G > £ex,* ~ 1 eV 
for yo = 3 and L ex > 10 9 L 0 /4”\ where R a = 7 x 10 10 cm and L G = 
4x 10 33 erg/s stand for the solar radius and luminosity, respectively. 
We discuss each of the requirements below. 

The typical seed photon energy lies in the optical energy range 
for mildly relativistic electrons, while the injection of relativistic 
electrons with y > 10 would push the typical energy below the 
far-IR regime. Such low energy photons are not, however, com¬ 
mon in GRB environments. For example, in the collapsar model 
for GRBs {Wooslevll 19931) . the burst is related to the core-collapse 
of a m assive star, typ i cally that of a Wolf- Rayet star (e . g. GR B 
980425 iGalama et all (1 19981) : GRB 030329 iHiorth et al.1 l2003h ). 
whose emission usually peaks in the far-UV/soft X-ray regime. A 
companion star would typically be of similar size and hence its 
emission is also expected to be dominated by this energy band. 

The fact that r ex » R a suggests that the seed photons for 
EC scattering cannot be directly provided by the progenitor star. 
A spatially extended source is, therefore, more plausible. One pos¬ 
sibility is that photons from the collapsar are reflected at the strong 
wind of the massive star. In this case, the energy density of the re¬ 
flected photons would still appear boosted by a factor of T 2 in the 
comoving frame of the emitting region. However, for typical val- 
u es fo r the mass loss rate and the wind velocity of the progenitor 
{Chevalier & Lill2000l) . the Thomson optical depth for scattering is 
given by 


_ cnM w ~ 3 Mw 10 8 cm/s 10 14 cm 

Tt 4/rr ex m p V w X 10- 5 M o /yr V w r ex 

which makes this scenario less efficient. Even if there were param¬ 
eters that led to tt ~ 1, the required luminosity of the external pho¬ 
ton field would still be many orders of magnitude larger than typi¬ 
cal luminosities of massive stars. The cocoon’s emission is stronger 
and even with such a small optical depth enough energy might be 
reflected back to serve as seed for external inverse Compton. How¬ 
ever, the typical energy of the cocoon’s photons is — 1 keV, which 
is much too high to serve as a seed for scattering to the sub-MeV 
y-rays. 

The luminosity constraint is relaxed, if we assume that 
4” J » 1. In principle, this is a viable assumption, but pushes the 
source into the Compton catastrophe limit. Roughly speaking, we 
showed that for » l/4y 0 , the luminosity of the (i + l)-th SSC 
generation becomes larger than preceding one (see Figs. HI and [3}. 
and eventually leads to L“£? N » L^S. Since, the peak energies of 
the two components differ, the final y-ray emission will be radically 
d iffer ent from the Band-like GRB observed spectrum {Band et alj 
Il993l) . More specifically, the electron compactness that leads to 
I* Bl > f-EC should be higher than 


r 


4yo 



2 r 4 L ex \ 1/JVT 

/ 


(65) 


emission is also challenged on a more general ground: the up-scattering of 
a black-body seed photon field by cooled electrons (n e oc y~~) results in 
much softer y -ray spectra than typically observed (e.g. iPreece etafll 19981 

l200d) . 



Figure 10. Photon spectra in 1 keV-1 GeV energy range obtained in the 
EC scenario for the prompt GRB emission. Different curves correspond 
to different values of the electron compactness (see legend). The arrow 
denotes the position of the observed peak energy, which is taken to be 
6 p obs = 0.2 MeV. For the rest of the parameters, see text. 


where we used eqs. {43}. {44} and {59}. The above relation may 
also be written as 



where we made use of eq. < 162 b and relation ( 163b . Even if Nj = 1, 
the numerical factor at the right hand side of the relation ( 166b is 
of order unity. Thus, if we try to relax the luminosity constrain by 
using £' e nj » 1, we find that the source is driven into the Compton 
catastrophe regime and at the same time LJ c s Nt » L^. 

We demonstrate the onset of the inverse Compton catastrophe 
through an indicative example, where we used the following param¬ 
eters: z = 1, Tk = 10 53 erg/s, eg = 10 -5 , T = 100, r em = 10 14 cm, 
4 — 4 m/r = 10 12 cm, and yo = 3. We considered a fiducial pho¬ 
ton field with te X ,* = 10 eV, r tx ~ r em , and L cx = 3 x 10 44 erg/s. 
The magnetic and external photon compactnesses for these values 
are £ B = 0.01 and 4x = 0.07. For the electron compactness, we 
adopted three values ranging from 10~ 2 up to 1; given the other pa¬ 
rameters, the maximum value of results in Lec ~ 10 49 erg/s (see 
eq. {44})- The high-energy spectra, which correspond to one GRB 
pulse with duration St ~ r b /cF ^ 0.3 s, are shown in Fig. [TO] For 
4” J <0.1, the EC component dominates the y-ray emission, at the 
cost, however, of low luminosity. Increase of 4" J transfers the y-ray 
power from the EC to the highest SSC component. Thus, although 
the power is ~ 10 51 erg/s for 4" J = 100, it is emitted at » MeV 
energies. Notice that we have already used a high fiducial value for 
f'ex- 

Concluding, we have shown, using simple arguments, that 
there are no reasonable physical parameters that lead to L ~ 
10 51 erg/s and at the same time suppress the Compton catastrophe 
effect. 
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5 DISCUSSION 

The inverse Compton catastrophe is essentially a runaway that 
transfers radiative power from low energy to high energy photons. 
It is applicable to leptonic plasmas where relativistic electrons ra¬ 
diate seed photons by synchrotron and/or inverse Compton scatter¬ 
ing, followed by further inverse Compton scatterings. If the source 
is compact enough in relativistic electrons, one can find parame¬ 
ter regimes, where the energy density of synchrotron photons be¬ 
comes larger than the magnetic one. This sets a threshold condition 
above which each successive generation of radiated photons due to 
the inverse Compton process becomes more important to electron 
losses than its previous one. This is a non-linear process, since it 
depends on internally produced photons, and could lead into an ex¬ 
ponentiation of electron losses, hence the term ‘inverse Compton 
catastrophe’. 

The aim of this paper is to examine the physical requirements 
for an optically thin to Thomson scatterings source to enter the in¬ 
verse Compton catastrophe regime and to calculate the expected 
radiative signatures in such a case. We have considered a one-zone 
model, where relativistic electrons are injected at a constant rate 
into a magnetized spherical source that is immersed inside an ex¬ 
ternal photon bath. We calculated, using standard expressions for 
the energy losses and radiative processes in the relativistic regime, 
the evolution of the system until it reaches a steady state. Because 
of the intrinsic non-linearity of problem, we restrict our analytical 
investigation to asymptotic regimes. The results of this analysis are, 
however, enlightening, and in rough agreement with our numerical 
ones, which are based on the simultaneous solution of the coupled 
kinetic equations for electrons and photons. This allows us to treat 
in detail the non-linear effects, which become central during the 
transition to the inverse Compton catastrophe regime. 

There are various source parameters that play an important 
role for the onset of inverse Compton catastrophe. As can be de¬ 
duced from the analytical treatment presented in Sect. 2, the key 
role belongs to the compactness of the injected electrons. An in¬ 
spection of eq. 0 reveals that as long as this is smaller than the 
magnetic compactness or the external photon one, then the sys¬ 
tem is linear and its solution leads to the standard expressions for 
the electron and photon distributions in the ‘fast’ or ‘slow’ cool¬ 
ing regimes. On the other hand, if the electron compactness in¬ 
crease, an increasing part of the losses is shifted to the internally 
produced synchrotron and inverse Compton photons. We have set 
the threshold for the inverse Compton catastrophe to be that value 
of the electron compactness that makes the electron losses due to 
external agents (magnetic and external photon fields) to be equal to 
the internal ones (synchrotron and Compton). Above that value the 
system becomes non-linear. 

A more sublime parameter for the onset of the inverse Comp¬ 
ton catastrophe is the characteristic injection energy of the electrons 
(yo), under the assumption that these are mono-energetic. This is 
related to the fact that the inverse Compton catastrophe enters its 
full modus operandi only when in the Thomson regime. The expo¬ 
nential photon growth in the source is naturally suppressed by the 
Klein-Nishina effects. Thus, as eq. 0 reveals, only small values 
of yo can produce several Compton generations before the Klein- 
Nishina cutoff is reached. 

As we have shown in Sect. 2, one can find parameters, albeit 
limited, that make the system enter the inverse Compton catastro¬ 
phe limit, while the electrons do not cool efficiently (slow cool¬ 
ing regime). This is in contrast to the standard notion, according to 
which, the inverse Compton catastrophe is alway associated with 


efficient electron cooling. So, at least according to our definition, 
the inverse Compton catastrophe is a misnomer because it does not 
necessarily imply catastrophic electron losses. 

An interesting feature of a source that is deep in the non-linear 
regime, i.e. when the electron injection compactness dominates 
over all the other compactnesses, is its multi-wavelength photon 
emission. In Sect. 3 we have shown that the photon spectra obtain 
an almost universal broken power-law shape, and it extend up to the 
maximum electron energy - see Fig.|7] The break occurs at ~ m e c 2 
in the rest frame of the source and is related to the Klein-Nishina 
suppression. The spectral index below the break energy does not 
depend on the power-law index of the injected electron distribu¬ 
tion, as is the case in the linear regime (both in the slow and fast 
cooling limits). In Sect. 2 and 3 we have shown that the spectral in¬ 
dex in the Compton catastrophe regime depends on the electron and 
magnetic compactness, but only logarithmically. Strictly speaking, 
the photon spectra in the Compton catastrophe regime do not satu¬ 
rate. Yet, they are largely independent of the specifics of the exter¬ 
nal photons, while their dependence on the electron and magnetic 
compactnesses is very weak. Thus, one can draw some analogies 
to the thermal Comptonization process, where the photon spectrum 
saturates into a specific power-law up to the maximum possible en¬ 
ergy. 

The multi-wavelength spectra from astrophysical sources, 
such as blazars and GRBs, do not generally resemble these ex¬ 
pected in the Compton catastrophe regime, which as we have 
shown bear some unique features. Assuming that the multi¬ 
wavelength emission of compact non-thermal sources originates 
from one region one can then apply the principles of the inverse 
Compton catastrophe to constrain their parameter space. The need 
of large luminosities when combined with the inferred small radii 
push the compactnesses of these sources to high values, and in or¬ 
der for the Compton catastrophe to be avoided, special conditions 
are required. We present two indicative applications in Sect. [4] For 
the case of blazars we set an upper limit to the size of the BLR re¬ 
gion. In the GRB case, we have shown how one can use the Comp¬ 
ton catastrophe to make a rather robust argument about ruling out 
an external Compton origin of the prompt (sub) MeV emission. 

We have studied the radiative signatures of compact sources 
in the inverse Compton catastrophe regime in a time-independent 
way by focusing on the steady-state photon spectra. These obtain 
a power-law shape that is the combined result of multiple inverse 
Compton scattered photon generations. As the flux of each under¬ 
lying photon generation depends on the various parameters in a 
different way, time variations in one or more parameters, such as 
the luminosity of injected electrons, will have a different impact 
on their respective emission. The end result of the escaping photon 
spectrum in such cases is non-trivial and requires a time-dependent 
approach that will be the subject of a future study. 
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APPENDIX A: ENERGY DENSITY OF HIGHER ORDER 
SSC PHOTON GENERATIONS 


We calculate the energy density of the i-th SSC photon generation 
(“ssc,i) under the assumption of inverse Compton scattering in the 
Thomson regime. A rough estimate of the number (At) of the SSC 
photon generations that will be produced in the Thomson regime 
can be found by 


3 

Tv-t < 4. 

where yi is defined as 


705-1 

m e c 2 



7o +1 > 


(Al) 


(A2) 


where we inserted ejm e c 2 = (4/3 )'byf' +l \ Thus, 


At = 



(A3) 


where the brackets denote the integer value of the enclosed expres¬ 
sion. Assuming B = 1 G and y 0 = 10 we find At = 5, whereas for 
more energetic electrons, e.g. yo = 10 3 , this reduces to At = 1. 

Under the (/-function approximation f or the synchrotron emis- 
sivity (see e.g. iMastichiadis & Knf]| 19951) . the number density of 
synchrotron photons is given by 


n s (x) = 


- l( B , -3/2 -1/2 

-r b cr T - -b x ' n t 

3 m e c 



(A4) 


where x is the photon energy in m e c 2 units. Using the same ap¬ 
proximation, the inverse Compton emissivity of a single electron is 
written as 


7ssc,i(-^ssc,i) — ^ (r T CM SSCi ,_i y 6 



(A5) 


and energy density of the i-th SSC photon generation is then given 
by 


^ssc,i — ^ 0"T^b 


J* dx S sc,i ^ 


dyy 2 n t {y)u ssc ; l -id\x i 


-y 2 x ssci -\ (A6) 


The stationary electron distribution n c is written as n e (y) = k s 6(y - 
yo) and k f y~ 2 for y < y 0 , in the limiting cases of slow and fast 
cooling, respectively. The energy density of synchrotron photons 
and successive SSC generations is then given by 


4 


^syn — 

j 

(4 \ 2 

^ssc,l — j 

[-o-T/yGej 

j 

(4 \ 3 

^ssc,2 = U B | 

[-O-TAGej 


(A7) 

(A8) 

(A9) 


u 


ssc,i 


U B | 2^TBj G t 


(A10) 


where 

G e = [ ks{ fO^ ^ SlOWC °° ling (All) 

| 4(4 J ,4x,4)yo, fast cooling 

Equations ( IA7MA101 are similar in both cooling regimes. All the 
information about the evolution of the electron distribution to a 
steady-state under the influence of escape, injection and cooling 
enters only through k s and k f , which are themselves functions of 
the basic compactnesses that describe the physical system. 


APPENDIX B: INCLUSION OF T HE EC ENERGY 
DENSITY IN ELECTRON COOLING 


The average energy of once inverse Compton scattered external 
photons, as measured in the rest frame of the emission region, is 
£ec = (4/3)yQ(Te CXi *), where yo is the electron Lorentz factor, T 
is the bulk Lorentz factor of the emission region, and e ex ,* is the 
characteristic energy of the external photons as measured in the 
rest frame of the galaxy. From this point on, we use unprimed and 
starred quantities, when these are measured in the rest frame of 
the emission region and the galaxy, respectively. In the above, we 
assumed that boosting of the external photon energy by a factor 
of ~ T because of the relativistic motion of the emission region 
with respect to the external source of photons. The inverse Comp¬ 
ton scattering of the EC photons by electrons with Lorentz factor 
yo = 100 takes place in the Klein-Nishina regime, since 


7o£ec ~ 10 7 eV 


/ yo y T 6 eXi * 

Vioo/ ToTeV 


» m e c 2 , 


(Bl) 


and one can safely neglect it from the electron cooling term (see 
eq. Oil). However, for lower Lorentz factors, e.g. yo = 10, the above 
relation becomes yofEC < meC 2 ', the first inverse Compton scat¬ 
tering of EC photons would occur in the Thomson regime, while 
the second scattering would marginally fall in the Klein-Nishina 
regime. We can therefore include only the first generation of EC 
photons in the electron cooling term. In what follows, we will per¬ 
form a similar analysis to the one presented in Sects. 2.2.1-2.2.2, 
but with the addition of the energy density of EC photons (mec) to 
the total one. 

Under the (5-function approximation for the inverse Compton 
emissivity (see Appendix[A}, the energy density of the EC photons 
is given by 


mec - -o’ji'bGcU^, (B2) 

where G e is defined in eq. dA lib . The similiraty to the energy den¬ 
sity of synchrotron photons is evident, with the external photon 
source playing, in this case, the role of the magnetic field (eq. lA7l l. 
Finally, one can define the respective photon compactness as 


4 

tEC - -O'T/yGetex- 


(B3) 


Bl Slow cooling regime 

The steady-state electron distribution is not determined by energy 
loss processes in this regime. Thus, the inclusion of uec in the elec¬ 
tron cooling term should not alter the solution for ;z e , which is writ¬ 
ten as n c (y) = k s 6(y - yo) with k s given by eq. ®. 

The total compactness of the source, in this case, is written as: 

Nt~ 1 

Act = 4 + 4x + 4(4 C j 7o) + 4x(4C J yo) + 4 (44 nj 7o) !+1 , (B4) 

i=l 












where we made use of eqs. ( 110b . ( II lb and ( IB3I) . Similar to what we 
have shown in Sect. 2.2.1, we may approximate the total compact¬ 
ness of the source as 


/ 4 + 4x, C<l/4yo 

“* ~ \ 4 (< nj ro) + 4 KVo) T , 4 nj > 1 /4yo, 

where for the latter, we assumed that 4 > 4 (44 70 ) ■ For 

4 J > l/4yo> the total compactness of the source is larger than the 
respective one in eq. 03- as expected. The slow cooling condition, 
i.e. £ > 1 or 3/(4/. nj yo) > 1, is satisfied if 


c < 

/ex < 

or 

c > 

/ex < 


1 

-and 

4yo 

3 


44ro ’ 



(B6) 

(B7) 


(B8) 

(B9) 


where / ex is defined in eq. {H- We note that the above constraints 
are identical to these of Sect. 2.2.1 apart from the last one (see 
eqs. d and {B9}. Inclusion of the uec makes the upper limit on 
/ ex more stringent than before. Thus, for the same 4 J and l B , the 
external photon source should be less compact, i.e. the ratio 4x/4 
should be smaller, in order to suppress electron cooling. 


B2 Fast cooling regime 


In the fast cooling regime, the steady-state equation of electrons 
may be written as 


k f 


N f—1 

/ ex (l +akt)+ (akf) i+l 

x 1=1 


9g j /4 
4o- T r b y 0 ’ 


(BIO) 


where a is defined in eq. {22}. In the above, we have also inserted 
the compactness term related to the EC photons, i.e. Iec = 4<*4- 
Same as in Sect. 2.2.2, eq. ( IB 1 01 1 can be studied in two regimes 
according to the value of akf. 


• for ak { < 1, the main energy loss channels for electrons are 
synchrotron radiation and inverse Compton scattering on external 
photons, synchrotron and EC photons. The total compactness can 
be approximated as 


Aot ~ 4 + 4x + Ay n + 4c - 4 /ex( i + akf). (Bll) 

The steady-state equation dB 101 then simplifies into a second order 
polynomial equation with the following solution: 


kf 


1 

2a 


-1 + 



v 

This is valid as long as akf < 1 or equivalently if 


(B12) 


/ex >6-^. (B13) 

t-B 

Equations dB 1 2b and dB13t should be compared to eqs. {23} and 
{26}. respectively. Taking also into account the fast cooling condi¬ 
tion, i.e. £ < 1, we find 


/ex > 


3 

44yod +4C j To) 


(B14) 
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Conditions d B 13b and dB 14b set a lower limit on / ex , which is writ¬ 
ten as 


/ ex > max 


44ro(i + 4CVo) 4 ) 


(B15) 


This is only slightly different than the one derived in Sect. 2.2.2 
(see eq. {33}). 

• for akf > 1, electrons may cool on the EC photons as well 
as on higher SSC photon generations. We are interested in the 
case where the 4 c,x t -i > 4c- In this asymptotic limit, the multi¬ 
wavelength photon spectrum will be dominated by the highest SSC 
photon generation. Thus, the features of the Compton catastrophe 
will not be hidden by a strong EC component. We may then write 


4t ~ 4 c,n t -i - 4 ( a k f ) Nj . 


(B16) 


In this limit, the energy density of EC scattered photons should 
not have an obvious effect on the electron steady-state distribution. 
We find, indeed, that the solution to the steady-state eq. ( IB 1 0I > is 
identical to eq. {28}, namely 



a 



1/(JV T +1) 


(B17) 


The conditions 4c, n t -1 > 4c or, equivalently, 4 < 4 (ak f ) NT 1 
and £ < 1 translate into the following two constraints on / ex and 

C s : 


/ex 


C 


< 


> 


1 + 


4 J 


Nj -1 
(Vy + 1 


— 

J 



(B 18 ) 

(B19) 


where the upper limit on / ex is slightly more stringent than before 
the inclusion of uec- 


Summarizing, we showed in detail that the inclusion of the energy 
density of EC scattered photons in the electron kinetic equation 
has an effect on the steady-state solutions derived only in the fast 
cooling regime and on the respective constraints. This is not un¬ 
expected, since electron cooling becomes important, by definition, 
in the fast cooling regime. Quantitatively speaking though, these 
effects are not important. 


APPENDIX C: THE COMPTON Y PARAMETER 

In Sect. [2] we derived the conditions, in terms of the various com¬ 
pactnesses, that result in the dominance of the inverse Compton 
emission over the synchrotron one. However, the strength of the in¬ 
verse Compton emission is traditionally expressed in terms of the 
Compton Y parameter. 

In the relativistic limit, the Compton Y parameter is defined as 
dRvbicki & Li ghtmaiill 19791) 

Y= ^(yo>max(T e ,T^), (Cl) 

where |(7 q) is the average energy gain per scattering and r c is the 
optical depth for Thomson scattering, which is given by 

T e = Cr T Wein,. (C2) 

In the above equation, n A = yon e (7o) an d is the only quantity enter¬ 
ing in the definition of Y that depends on f"' 1 and C B . In the slow and 
fast cooling regimes, we find, respectively, that n e i = k s and kf/y q. 
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Here, k s and k f are given by eqs. ®, (El and ( 128b . respectively. 
The Thomson optical depth can be, thus, written as 

—— slow cooling 


|-1 + -^1 + -^r^j fast cooling 

fast cooling (CC limit) 


(C3) 


If r e » 1 (or, in practice, r c > 10) our analysis is not valid, 
since other physical processes that we do not take into account, 
such as photon-photon absorption, become important. It is inter¬ 
esting though, that even for r e < 1-3, there are parameters that 
drive the system to inverse Compton catastrophe. For simplicity, 
let us consider the slow cooling case. In this case, higher order SSC 
photon generations become more luminous for > l/4yo (see 
Sect.d, which translates into r c > 3/4. Thus, in the slow cool¬ 
ing regime and for 3/4 < r e < 1 - 3 the system may be driven to 
the Compton catastrophe regime, while being still optically thin to 
scatterings. 

For the parameter regime that does not lead to the inverse 
Compton catastrophe and for r e < 1, the Y parameter is written 
as 


44 nJ To sl° w cooling 

^ |-1 + yj 1 + fast cooling 


(C4) 


Note that both expressions can be derived by Y = £i C /£ S yn ~ 
Asc,i/4 y n (see eqs. ( lilt . 110b . 05b and P4b ): the approximation 
breaks down in the Compton catastrophe regime. In the absence of 
external radiation fiels (/ cx = 1), the Y parameter in the fast cooling 
regime reduces to 




12t“ J 

— 

124 nj 


» 1 


<K 1 


(C5) 


By identifying 12^ as e e and t B as e%, the above expressions 
are equivalent to these commonl y used in GRB literature (e.g. 
ISari et alJl 19961 : ISari & Esli]|200ll) when the scatterings take place 
in the Thomson regime (f or the expression of Y in the presence of 
Klein-Nishina effects, see lNakar et alj|200lj) . 


APPENDIX D: OBSERVED EC FLUX FOR 
MONO-ENERGETIC ELECTRONS 

In the slow cooling regime, the electron distribution is defined by 
the balance between the injection and escape terms. As long as the 
escape term is independent from energy, the stationary electron dis¬ 
tribution has the same energy dependence as the injection term. In 
our analysis we focus on mono-energetic injection. 

Assuming that the electron distribution in the comoving frame 
is also isotropic we may write n c (y, £1') = k s /(4n)6(y - yo), where 
primed quantities are measured in the blob’s frame, while quantities 
with the subscript “obs” are measured in the observer’s frame. By 
inserting the above expression into eq. (5) in D95 we find the EC 


emissivity in the comoving frame: 

■ , / rvy CCTMex^To / ^ r «ex,*(! + H') y - ,\ m n 

Jec (e , £2 ) = — | j |---) 8(e - e 0 ) , (Dl) 

where we used the property d(g(x)) = 2,d(x - xi)l\g\xi)\ with x,- 
being the roots of g{x) = 0, and 

fo = yoTeex,*(l (D2) 

Substitution of )ec (e', £2') into eq. (1) in D95 results in 

r’obs/ r-> x _ r- / 2)( 1 + rt,bs) j / / e 0 bs(l + z) \ / r x 

(Cobs; Hobs) — Fq I -j - —-I I -I 8 (fobs ~ Cec) • (D3) 

where we made use of the relation 

p/ J + ^obs - P \ _ 1 + /Jobs 

\ 1 - P/J-obs ) 1 + P 

and 

j-, D 4 V b o- T c 

To = ————u ex K s yo 


AnDl 


£ec = 


2? 2 e e x,*ro(l +/iob S ) 

(1 +z)(l +P) 


(D4) 

(D5) 

(D6) 


APPENDIX E: EXPRESSIONS FOR THE LOSS AND 
INJECTION OPERATORS USED IN THE NUMERICAL 
CALCULATIONS 

The numerical results presented in Sect. [3] were obtained using a 
numerical code that solves the coupled partial differential equations 
that govern the evolut ion of photons and relativi stic electron^/] in 
both energy and time dMastichiadis & Kirkll 19951) . It is more con¬ 
venient to rewrite eqs. ( 153b and ( 154b in dimensionless form, where 
time is in units of the crossing time (r = tc/r b ), the photon and elec¬ 
tron energy is in units of m c c 2 and the particle number densities are 
normalized as 

«ei(y,T) = nPy,T)cr T r b (El) 

hy(x, t) = n 7 (x,T)cr T r b , (E2) 

where <xx is the Thomson cross section. In what follows, we 
present the expressions of the loss and source operators that appear 
in eqs. El-Eli, after transforming them into a dimensionless 
form. 

• Synchrotron radiation 


QT(x,t) = 

J' dy h c (y,T)J syn (x,y), 

(E3) 

I7 n (y,r) = 

^ I«e (y. T) J' dx] syn (x, y)j , 

(E4) 

where the single particle synchrotron emissivity is 


7*syn(-^»y) — 

V3 q\B sin 9 r b / x \ 
hxm c c 2 c \x c ) 

(E5) 

Xc = 

3 . 2 

— sin 6by~, 

(E6) 


where F(z) is defined in eq. (6.31c) of iRvbicki & Lightmanl 
H979t) . 6 is the pitch angle between the magnetic field line and the 


9 We do not distinguish between electrons and positrons. We treat them as 
one particle population and we commonly refer to both as electrons. 
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particle’s velocity, and b = B/B ct . 
• Synchrotron self-absorption 


We treat synchrotron self-absorption only as a sink term for 
photons, while we neglect the respective energy gain term in the 
electron equation. The loss term in the photon equation is written 
as 


_£“ a O, r) = Ra ssa h r {x, r). 


(E7) 


In the above equation o- ss a is the synchrotron self-abso rption coeffi¬ 
cient (c.f. eq. (6.50) from lRvbicki & Lightmanl il979l) ) 


^ssa 


C 3 h 3 1 r° 2 _fl_ I h c (y, t) 
8^(m e c 2 ) 3 ct t R 2 x 2 77 fly \ y 2 


%.yn(x, y) 


(E8) 


where j syn is defined in eq. iE5b . 

•Inverse Compton scattering 

The energy injection term for photons is given by 


q; cs (x„t) 

where 


Jflyn e (y,r) J dxj ics {x, jci, r), 


j\c S {x, xi ,r) = 


3 xi h y {x) 

4J — 


F c (q,G). 


(E9) 


(E10) 


In the above, function F ( is given by eq. (2.48) in 
iBlumenthal & Gould dl970bl) : 


Fc(q,G) = 2q\nq + (1 + 2q){\ - q) H— 


1 (GqT 


2 (1 + Gq) 


(1 ~q), (Ell) 


where 


G = 
q = 

Ei = 



G(l-£,) 

x\ 


7 


(El 2) 

(El 3) 

(E14) 


The electron energy loss term is split into two parts, = -£e CS " T + 
_£‘ cs,kn . For collisions taking place in the Thomson regime the loss 
term is given by 


“■ T (y,r) = “ (y 2 n e (y,T)U ph j) 

(El 5) 

r 3/4y 

U P hT (t ) = I dxxh y (x,T) 

Jo 

(El 6) 


while for scatterings occuring in the Klein-Nishina regime we may 
write 


f-.KN, T ) » 3 ,%(y ’ T) 

c V7 ’ ; 8 y 


J dx 

3/4y 


h y (x, r) 


ln(4xy) - y 


(E17) 


where we made use of eq. (2.57) in lBlumenthal & Gould d 1 970bl) . 
The above expression for the inverse Compton losses is based 
on the assumption that the electron losses all its energy in one 
collision while in the Klein-Nishina regime. 


•Photon-photon pair production 
The sink term for photons is written as 


-Cy y (x, t) = h y ( x, r) 


J dx' hy(x'. 

o 


T)Ryy(XX'), 


(El 8) 


where for the reaction rate R vv we use the approximate expression 
(4.7) of lCoppi & Blandfordl dl990h 

w 1 — 1 

Ryy(w)~ 0.652-— lnwtf(w-l), (E19) 

w 3 

where w = xx' and H(z ) is the usual Heaviside function. Assum¬ 
ing that the energy of the absorbed photon is shared between the 
electon and positron, the respective source term for pairs becomes 

= 4n 7 (2y, r) f dx' h y (x\ T)R ry (2yx'), (E20) 

Jo 


APPENDIX F: REPEATED SCATTERINGS OF A FIXED 
BLACK-BODY PHOTON FIELD 

We showed that, in the pure SSC case, the photon spectrum in 
the inverse Compton catastrophe limit can be described by a bro¬ 
ken power-law, with the spectral index below the break (/? cc ) de¬ 
pending on T’e" 3 , among other parameters. We found that typically 
ftcc ^ 0-6> while only for unphysically high £™ J , the photon spec- 
tra become harder. This is in contrast with one of the findings in 
IZdziarski & Lambl d 19861 ) - hereafter, ZL86, where photon spectra 
with/? ~ 0.3 were obtained for values lower than these adopted 
in this study. 

To test if there is a physical reason for this discrepancy we 
calculate the photon spectra using the same parameter values as in 
Fig. 1 in ZL86: B = 0, R = 8 x 10 6 cm, y min = 1, y mllx = 600, 
s = 2.2. Our definition for the electron compactness has an extra 
factor An in the denominator (see eq. 0 ). Thus, f'T = (/{An) = 2.4, 
for ( = 30. Moreover, the seed photon field is assumed to be a 
black-body with 7 bb — x SO ft/2.7 k =s 4 x 10 6 K and compactness 
(soft = {(Lsott/L) = 0.35, where L soft /L = 1/85. 

The comparison between the spectra is shown in Fig. ED The 
result of ZL86 is plotted with a black line, while our numerical 
results are shown with red and blue lines. The red curves are ob¬ 
tained usin g a simplified expression for the Compton scattering rate 
(eq. (44) in lMastichiadis & Kirkll 19951 ). which is more close to the 
one used by ZL86. We also neglected scatterings taking place in the 
Klein-Nishina regime. The blue curves are obtained using the full 
expression (in the relativistic limit) for the Compton scattering rate 
(see eq. ( IE9I ) ). Finally, to demonstrate the effect of photon-photon 
(yy) absorption on the spectra, we artificially deleted the respective 
terms from the electron-photon equations (dashed lines). A few re¬ 
marks follow: 

• for x < 1, all spectra have similar spectral index (~ 0.3). 

• the discontinuity found at x ~ 1 by ZL86 is a numerical arti¬ 
fact, that does not appear in our results. 

• for x > 1, our spectra (with yy absorption included) are 
steeper than the one in ZL86. We note that ZL86 used a more crude 
expression for the yy cross section than ours (see eq. <El9l ). 

• a comparison of the dashed lines and the black solid line 
shows that for x > 1 the red curves are more close to the result 
by ZL86, since they are obtained using a similar expression for the 
Compton scattering rate. 

Concluding, there are some differences at the high-energy part of 
the spectrum (x > 1), which can be understood in terms of differ¬ 
ences in the yy cross section and in the Compton scattering rate. 
Yet, there is good agreement at the low-energy part of the spectrum 
(x < 1). This is harder than the spectrum we typically obtain in the 
Compton catastrophe regime, whenever the synchrotron photons 
serve as the seed photons for scattering. 
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log X 


Figure FI. Photon spectra obtained for the same parameters used in Fig. 1 
by ZL86. The spectrum of ZL86 is shown with a black line, while our results 
are shown as red and blue curves. Different lines have following meaning: 
simplified expressi on for the IC scattering rate (r ed lines), exact expression 
for the IC rate from lBlumenthal & GoulJ1970bl (blue lines), yy absorption 
included (solid lines), yy absorption neglected (dashed lines). 
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